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We introduce and analyze a class of one dimensional insulating Hamiltonians which, when adi- 
abatically varied in an appropriate closed cycle, define a "Z2 pump". For an isolated system a 
single closed cycle of the pump changes the expectation value of the spin at each end even when 
spin orbit interactions violate the conservation of spin. A second cycle, however returns the sys- 
tem to its original state. When coupled to leads, we show that the Z2 pump functions as a spin 
pump in a sense which we define, and transmits a finite, though non quantized spin in each cycle. 
We show that the Z2 pump is characterized by a Z2 topological invariant that is analogous to the 
Chern invariant that characterizes a topological charge pump. The Z2 pump is closely related to 
the quantum spin Hall effect, which is characterized by a related Z2 invariant. This work presents 
an alternative formulation which clarifies both the physical and mathematical meaning of that in- 
variant. A crucial role is played by time reversal symmetry, and we introduce the concept of the 
time reversal polarization, which characterizes time reversal invariant Hamiltonians and signals the 
presence or absence of Kramers degenerate end states. For non interacting electrons we derive a 
formula for the time reversal polarization which is analogous to the Berry's phase formulation of the 
charge polarization. For interacting electrons, we show that abelian bosonization provides a simple 
formulation of the time reversal polarization. We discuss implications for the quantum spin Hall 
effect, and argue in particular that the Z2 classification of the quantum spin Hall effect is valid in 
the presence of electron electron interactions. 
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I. INTRODUCTION 

In recent years, the advent of spintronics has motivated 
the search for methods of generating spin currents with 
little or no dissipation. One class of proposals involves 
designing an adiabatic pump in which the cyclic variation 
of some control parameters results in the transfer of spin 
across an otherwise insulating structurei*2iM. Such a 
spin pump has been realized in quantum dot structures^. 
A second class of proposals involves using the spin Hall 
effect to generate a spin current using an electric field®-^. 
Interest in this approach has been stimulated by the ex- 
perimental observation of spin accumulation induced by 
the spin Hall effect in doped GaAs structures®'^. In these 
experiments the spin current is accompanied by a dissipa- 
tive charge current. This motivated Murakami, Nagaosa 
and ZhangiS to propose an interesting class "spin Hall in- 
sulator" materials which are band insulators that have, 
according to a Kubo formula, a large spin Hall conductiv- 
ity. However, the spin current which flows in the bulk of 
these materials is not a transport current, and can not be 
simply measured or extracted. A crucial ingredient for 
the generation of transport currents is the existence of 
gapless extended edge states. Such states are generically 
not present in the spin Hall insulatorsii. 

Motivated by the spin Hall insulator proposal, we in- 
troduced a model of graphene in which the symmetry 
allowed spin orbit interactions lead to a quantum spin 
Hall effectiSii^. A related phase has been proposed for 
GaAs in the presence of a uniform strain gradienliM. This 
phase is characterized by a bulk excitation gap and gap- 
less edge excitations. In the special case where the spin 
Sz is conserved, this phase can be viewed as two copies 



of the quantum Hall state introduced by Haldanei^. The 
phase persists, however, in the presence of spin noncon- 
serving interactions as well as disorderiiSiiii^. Time re- 
versal symmetry protects the gapless edge states when 
electron interactions are weak, though strong interactions 
can open an energy gap at the edge accompanied by time 
reversal symmetry breaking^^'^^-. We argued that the 
quantum spin Hall phase is distinguished from a band 
insulator by a Z2 topological inde»13,, which is a prop- 
erty of the bulk system defined on a torus. We suggested 
a formula for this index in terms of the Bloch wavefunc- 
tions. However, the physical meaning of this formula and 
its relation to the edge states was not explicit. 

When placed on a cylinder (or equivalently a Corbino 
disk) , the quantum spin Hall system defines a kind of adi- 
abatic pump as a function of the magnetic flux threading 
the cylinder. In the case where Sz is conserved, advanc- 
ing the flux by one flux quantum results in the transfer of 
spin ASz — h from one end of the cylinder to the other. 
This is a spin pump, whose operation is analogous to a 
charge pump which could be constructed with a quantum 
Hall state on a cylinder. As envisioned by Thouless and 
CO workers in the 1980'3i2i2S, the adiabatic charge pump- 
ing process is characterized by a topological invariant - 
the Chern number - which is an integer that determines 
the quantized charge that is pumped in the course of a 
cycle. Equivalently, the Chern number provides a topo- 
logical classification of the two dimensional quantum Hall 
g-(-g^^g2ij22j23j24^ When Sz is conserved, similar ideas can 
be used to describe a quantized adiabatic spin pump*. 

A local conservation law is essential for this type of 
topological pumping process. For a finite system with 
closed ends, the eigenstates before and after a complete 
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cycle must be distinct. This means that two energy levels 
must cross in the course of the cycle. In the case of the 
charge pump, that level crossing is protected by local 
charge conservation because the two states differ in their 
eigenvalue of the charge at each end. In the absence 
of a conservation law there will, in general, be no level 
crossings, and the system will be in the same state before 
and after the cycle. 

Unlike charge, spin does not obey a fundamental con- 
servation law, so unless spin non conserving processes can 
be made very small it is not possible to define an adia- 
batic spin pump that works analogously to the Thouless 
charge pump. Nonetheless, in Ref. [T^ we argued that 
time reversal symmetry introduces a conservation law 
which allows for a topological pumping process. Specif- 
ically, we showed that for a quantum spin Hall state on 
a cylinder the eigenstates before and after adiabatic flux 
insertion are orthogonal, and can not be connected by 
any local time reversal invariant operator. When a sec- 
ond flux is added, however, the system returns to its 
original ground state. In this sense, the quantum spin 
Hall effect defines a pump" . In one cycle there is no 
charge transferred between the two ends. Since spin is 
not conserved, the two states can not be distinguished by 
a spin quantum number (though the expectation value of 
the spin at the end changes by a non quantized amount). 
The question therefore arises: What is it that is pumped 
between the ends of the cylinder? 

In this paper we examine this issue carefully and in- 
troduce a class of one dimensional models that exhibit a 
similar pumping behavior that is protected by time re- 
versal symmetry. The "Z2 spin pump" is analogous to 
the quantum spin Hall effect in the same sense that the 
charge pump is analogous to the ordinary quantum Hall 
effect. We introduce the concept of the time reversal po- 
larization, a Z2 quantity which signals whether or not 
a time reversal invariant one dimensional system has a 
Kramers degeneracy associated with its ends. We show 
that the change in the time reversal polarization in the 
course of an adiabatic cycle is related to a Z2 topologi- 
cal invariant which distinguishes a Z2 spin pump from a 
trivial cycle of an insulator which pumps nothing. This 
Z2 invariant is equivalent to the invariant introduced in 
Ref. ^^to characterize the quantum spin Hall effect. The 
present work, however, provides an alternative formula- 
tion which clarifies both the physical and mathematical 
meaning of the invariant. 

We study a family of one dimensional Hamiltonians 
which have a bulk energy gap and a length which is much 
larger than the exponential attenuation length associated 
with that gap. We suppose the Hamiltonian depends 
continuously on a "pumping parameter" t, satisfying the 
following properties: 

H[t + T]^H[t], (1.1) 

H[-t]^eH[t]Q-\ (1.2) 
where Q is the time reversal operator. In the case that 



the one dimensional system corresponds to a two dimen- 
sional system on a cylinder t/T may be viewed as the 
magnetic flux threading the cylinder in units of the flux 
quantum. In the course of the cycle time reversal sym- 
metry is broken. However, the second constraint ensures 
that the system passes through two distinct points t^ ^ 
and = T/2 at which the Hamiltonian is time reversal 
invariant. Condition H1.2|l may be relaxed somewhat, but 
it is essential that there exist two distinct time reversal 
invariant points tl and where (|1.2|l is locally valid. 
The existence of two such points plays a crucial role in 
the topological classification of the pumping cycle. In 
particular, we will show that pumping cycles in which 
H[tl] and H[t2] have different time reversal polarization 
are topologically distinct from trivial cycles. 

We will begin in section II by introducing a simple 
one dimensional tight binding model which exemplifies 
the Z2 spin pump. This model is closely related to a 
model of a spin pump that was recently introduced by 
Shindoui which may be applicable to certain spin 1/2 
quantum spin chains, such as Cu-benzonate and Yb4As3. 
This tight binding model incorporates spin non conserv- 
ing spin orbit interactions and provides a concrete illus- 
tration of the Z2 pumping effect. 

In section HI we provide a general formulation of the 
time reversal polarization for non interacting electrons. 
Our discussion closely parallels the theory of charge 
polarizationSSiSSiSLSSiSS, in which the charge polarization 
is related to the Berry's phase of Bloch wavefunctions. 
We show how the change in the time reversal polariza- 
tion defines a Z2 topological invariant characterizing the 
pumping cycle. 

In section IV we argue that the notion of time reversal 
polarization and the topological classification that fol- 
lows from it can be generalized to interacting systems. 
We describe an interacting version of our Id model us- 
ing abelian bosonization. This provides a different for- 
mulation of the time reversal polarization, which is well 
defined in the presence of interactions. 

In section V we conclude by addressing two issues. In 
VA we discuss the implications of the time reversal polar- 
ization for the quantum spin Hall effect. We argue that 
the two dimensional quantum spin Hall phase is a distinct 
phase from a band insulator even in the presence of elec- 
tron electron interactions. We then prove that this phase 
either has gapless edge excitations or exhibits a ground 
state degeneracy associated with time reversal symmetry 
breaking at the edge. We also comment on a proposal by 
Sheng, Weng, Sheng and Haldane'^° to classify the quan- 
tum spin Hall effect in terms of a Chern number matrix. 

In VB we ask whether the Z2 spin pump we have de- 
fined can actually pump spin. Despite the fact that the 
isolated pump returns to its original state after two cy- 
cles, we argue that when connected weakly to leads, the 
Z2 pump does pump spin, although the amount of spin 
pumped in each cycle is not quantized. 

In the Appendix we relate different mathematical for- 
mulations of the Z2 topological invariant. We begin by 
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showing that, hke the Chern invariant, the Zi invariant 
can be interpreted as an obstruction to globaUy defin- 
ing wavefunctions, 'provided a constraint relating time re- 
versed wavefunctions is enforced. We then prove that 
the invariant derived in this paper is equivalent to the 
one introduced in Rcf. 13. 



II. TIGHT BINDING MODEL 

In this section we introduce a one dimensional tight 
binding model of the Zi spin pump. This model is closely 
related to a model introduced by Shindou as an adiabatic 
spin pumpii. Shindou considered an antiferromagnetic 
spin 1/2 Heisenberg chain to which two perturbations 
which open a gap in the excitation spectrum are added. 
The first term is a staggered magnetic field hst which 
locks the spins into a Neel ordered state. The second is a 
staggered component to the exchange interaction AJst, 
which leads to a dimerizcd state. Interestingly, Shin- 
dou suggested that this model may be relevant to certain 
S" = 1/2 quantum spin chains, such as Cu-benzonate and 
Yb4As3, in which spins reside at two crystallographically 
inequivalent sites. He argued that in these systems h^t 
can be controlled by applying a uniform magnetic field, 
and AJst can be controlled with a uniform electric field. 

Shindou showed that a cycle in which AJ^t and h^t 
are adiabatically varied defines a topological spin pump, 
which transfers Sz = h in each cycle. The topolog- 
ical quantization of Shindou's pump requires the con- 
servation of Sz- In general, however, Sz non conserv- 
ing processes are allowed by symmetry. In particular 
the Dzyaloshinskii-Moriya interaction, d • (Si x S2), is 
allowed, and will inevitably lead to the violation of Sz 
conservation. We will argue, however, that provided the 
system retains time reversal invariance when hst = this 
system remains a Z2 spin pump even in the presence of 
the Dzyaloshinskii-Moriya interaction. 

We study a non interacting electron version of the 
Shindou model, where in addition to the spin degree of 
freedom we allow charge fluctuations. Consider a one 
dimensional tight binding model with a staggered mag- 
netic field, a staggered bond modulation as well as a time 
reversal invariant spin orbit interaction, 

H = Ho + Vh + Vt + Vso, (2.1) 

where 

Ho=toY^ (c^Q+la + 4+laC»a) , (2.2) 
i,a. 

yh = hstY.{-iyolpclc.p, (2.3) 

Vt = Atst 5](-l)XcLc«+ia + 4+i„C,„) (2.4) 




FIG. 1: (a) Schematic representation of the groundstate of 
Eqs. 12.112.51 for various t. The groundstates at the time 
reversal invariant points t = Q and t = T/2 are distinguished 
by the presence of Kramers degenerate end states, (b) Single 
particle energy levels E„{t) for a 24 site chain with Atgt/to = 
.4, h^lt/to ~ .8 and eso/to = .Ij/. (c) Single particle energy 
levels En{t) for a 24 site chain with 12 extra sites added at 
each end. (d) Low energy many body energy levels associated 
with one end of the chain. The degeneracy at t — T/2 and 
t = 3T/2 is protected by time reversal symmetry. The inset 
shows how electron electron interactions lifts the degeneracy 
at t = 0, f = r and t = 2T. 

and 

Vso = ■ ^af}{claC,+i0 - cl_^_^^Cif3). (2.5) 

Here Cgo is an arbitrary vector characterizing the spin or- 
bit interaction. This term explicitly violates the conser- 
vation of Sz, playing a role similar to the Dzyaloshinskii- 
Moriya interaction in Shindou's model. We consider an 
adiabatic cycle in which 

(At,t, /i,t) = (AtO, cos(2^t/T), /lO, sin(27ri/r)). (2.6) 

Since Vh is odd under time reversal, while Vt is even, 
condition (|1.2|l is clearly satisfied. At i = and T/2 the 
Hamiltonian is time reversal invariant. 

In Fig. 1(a) we depict groundstates in the strong cou- 
pling limit at representative points along the cycle. At 
t = T/4 and t = Vh dominates and locks the spins 

into a Neel ordered state. At t = and t = T/2, Vt dom- 
inates, and the system is dimerized with singlet pairs 
of electrons occupying alternate bonds. Importantly, the 
groundstate at i = T/2 is distinguished from the ground- 
state at i = by the presence of unpaired spins at each 
end. 

When Vso — Sz is conserved, and this model de- 
scribes a spin pump. In this case Vh + Vt can be de- 
composed into two independent periodic potentials which 
lock the densities of the up and down spin particles and 
slide in opposite directions as a function of t. As t evolves 
from to T, the periodic potentials slide by one lattice 
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constant. Provided there is "space" for the added spin 
at the ends, spin h will accumulate at the end following 
each cycle. 

We wish to understand how this spin pump is modified 
when Vso ^ 0, so that Sz is not conserved. In Fig. 1(b) 
we plot the single particle energy levels for a 24 site chain 
as a function of t for nonzero Vso- The bulk energy gap 
can be clearly seen with continuum states above and be- 
low. The energy levels that cross the gap are end states. 
Each line consists of two states which are localized at 
opposite ends. The crossing of the end states at i = T/2 
will play a critical role in what follows. When Vso — 
the degeneracy at T/2 is protected by spin conservation 
because the two states at each edge have Sz = ±h/2. 
Nonzero Vso does not lift the degeneracy provided the 
Hamiltonian remains time reversal invariant at T/2. The 
two end states form a Kramers doublet whose degeneracy 
can not be broken by any time reversal invariant pertur- 
bation. 

Because of the level crossing at t — T/2 it is clear 
that a system which starts in the ground state at f = 
will be in an excited state ai t = T. However, since 
the end states merge with the continuum the excitation 
will not be localized near the edge, and bulk particle- 
hole pairs will be excited. This is because there is no 
"space" to put the excitations at the ends. In section 
VB we will discuss the effect of connecting this pump 
to reservoirs which allow the end states to be "emptied" 
without exciting bulk particle- hole pairs. For the purpose 
of this section, however, we will study the operation of 
an isolated pump by adding several sites at the ends of 
the chain for which Vh and Vt vanish. This introduces 
additional midgap states localized at each end, allowing 
the cycle to proceed without generating bulk excitations. 

Fig. 1(c) shows the energy levels as a function of t with 
the extra sites added at each end. There are now several 
midgap states at each end. Since all of the midgap states 
are localized at one end or the other, the low energy 
excitations of the system can be factorized as a product 
of excitations at each of the two ends. 

In Fig. 1(d) we plot for < i < 2T the energies of 
the lowest few many body eigenstates associated with a 
single end, obtained by considering particle hole excita- 
tions built from the single particle states localized at that 
end. Though this picture was computed for non interact- 
ing electrons, it is clear that the Kramers degeneracy of 
the groundstate at t = T/2 and 3T/2 will be robust to 
the addition of electron electron interactions. The first 
excited state at T = 0, T and 2T in Fig. 1(d) is four 
fold degenerate (the middle level coming into that point 
is doubly degenerate). This degeneracy, however, is an 
artifact of non interacting electrons. The degeneracy is 
present because there are four ways of making particle 
hole excitations with two pairs of Kramers degenerate 
states. Electron electron interactions, however, will in 
general split this degeneracy, as shown in the inset, so 
there will be no level crossing at t — T. 

We thus conclude that when the isolated pump starts 



in its ground state at i = 0, it arrives in an excited state 
after one complete cycle at t — T. After a second cy- 
cle, however, at t = 2T the system returns to its orig- 
inal state. For this reason, we call it a "Z2 pump" . It 
is possible that by coupling to other degrees of freedom 
an inelastic process (such as emitting a phonon) could 
cause the excited state to relax back to the ground state. 
Nonetheless, there is an important distinction between 
this adiabatic process which generates an excited state 
and one that does not. In section VB we will return to 
this issue when we discuss connecting the pump to leads. 
The nontrivial operation of a single cycle depends criti- 
cally on the time reversal symmetry att = T/2. Breaking 
time reversal symmetry at that point leads to an avoided 
crossing of the energy levels, so that the system returns 
adiabatically to its original state at t = T. 

From the point of view of the end states, the non trivial 
pumping effect arises because there exist Kramers degen- 
erate end states at T = T/2, but not at T = 0. In the 
next section we show that this property is determined by 
the topological structure of the bulk Hamiltonian, H(t). 



III. TIME REVERSAL POLARIZATION AND Z2 
INVARIANT 

In this section we introduce the time reversal po- 
larization for non interacting electrons and show that 
changes in it define a topological invariant. Our dis- 
cussion will parallel the theory of charge polarization in 
insulator32Si2SiSL22i2S. In order to establish this connec- 
tion and to define our notation we will therefore begin 
by reviewing that theory, which relates the charge polar- 
ization to the average center of Wannier orbitals, which 
in turn are related to the Berry's phase of the Bloch 
wavef unctions. We next consider the role of Kramers' 
degeneracy in time reversal invariant systems and define 
a corresponding time reversal polarization in terms of 
the difference between the Wannier centers of Kramers 
degenerate bands. Finally, we show that the change in 
the time reversal polarization between t — and t = T/2 
of the pumping cycle defines a Z2 topological invariant 
which distinguishes a nontrivial Z2 pump from a trivial 
cycle. 



A. Review of theory of charge polarization 

Consider a one dimensional system with lattice con- 
stant a = 1, length L = Nc with periodic boundary con- 
ditions and 2N occupied bands. The normalized eigen- 
states for the n'th band can then be written in terms of 
cell periodic Bloch functions as 

IVn.fc) = ^KT^^'^M- (3.1) 
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We may define Wannier functions associated with eacfi 
unit cell associated with lattice vector R as 



\R, n) = — / dke 



-ik{R-r) 



Uk,n)- 



(3.2) 



The Wannier functions are not unique because they de- 
pend on a gauge choice for \uk,n) ■ In addition to changing 
the phases of the individual wavefunctions, the wavefunc- 
tions can be mixed by a general U (2N) transformation 
of the form 



\Uk.r. 



{k)\uk. 



(3.3) 



After this transformation, \uk,n) need no longer be the 
individual eigenstates of the Hamiltonian, but rather 
should be interpreted as basis vectors spanning the space 
spanned by the 2N occupied eigenstates. The Slater de- 
terminant of the 2N wavefunctions is unchanged up to a 
phase. 

Marzari and VanderbiltSS have provided a prescription 
for choosing Unm{k) to optimally localize the Wannier 
wavefunctions. Here, however, we are concerned with the 
total charge polarization, which is insensitive to the de- 
tails of Unm{k) ■ The polarization is given by the sum over 
all of the bands of the center of charge of the Wannier 
states associated with i? = 0, and may be writte n^^i^^ 



where the U{1) Berry's connection is given by 
A{k) 



Uk,n[^k\uk,n)- 



(3.4) 



(3.5) 



The integral is over the Brillouin zone from k — — tt to 
TT. If we require that the wavefunction |V'n,fc) be de- 
fined continuously in the reduced zone scheme, so that 
|V'n,-7r) = |^Tt,7r) then ^(— tt) = A{-k), and the integral 
may be considered to be on a closed loop, despite the 
fact that \un,k) is discontinuous from — tt to tt— . Under 
U (27V) a transformation which preserves this continuity 
Pp is invariant up to a lattice constant. For a transforma- 
tion in which the U (1) phase of Umnik) advances by 27rrn 
when k advances around the Brillouin zone Pp — > Pp + m. 
This reflects the fact that the polarization can only be 
defined up to a lattice vector. 

Changes in the polarization induced by a continuous 
change in the Hamiltonian H[t\ are, however, well de- 
fined. Thus, if the wave functions jw^.n (t)) are defined 
continuously between ti and ^2 for all k in the Brillouin 
zone, then we may write 



Ztt 



dk A{t, k) 



dk A{t,k) 



(3.6) 

*1(2)- 

Using Stokes theorem, this can be written as an integral 



where Ci(2) is the loop fc = — tt to tt for fixed t 
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FIG. 2: Schematic one dimensional band structure with spin 
orbit interactions. The energy bands come in time reversed 
pairs which are degenerate at = and k = n. 



of the Berry curvature 

.F(t,fc)=*^((Vt Uk,n{t)\'VkUk,n{t)) - C.c) (3.7) 
n 

over the surface T12 of the cylinder spanned by k and t 
bounded by ci and C2: 



Pph] - Ppih] = ^ I dtdk T{t, 



k). 



(3.8) 



For a periodic cycle H[t + T] — i?[t], the change in 
the polarization over one cycle, Pp{T) — Pp{0) is given 
by the integral in (|3.8|) over the entire torus defined by 
t and k. This quantity is an integer and defines the first 
Chern number associated with the wavefunction \uk,n{t)) 
on the torus. The Chern number characterizes the charge 
pumped in each cycle. For a cycle which satisfies the time 
reversal constraint in Eq. H1.2|l . T{—t,—k) = —J^{t,k), 
so the Chern number is equal to zero. 



B. Time reversal Polarization for Kramers 
degenerate bands 

Consider now a time reversal invariant system. The 
time reversal operator has the form 



(3.9) 



where Sy is spin operator and K is complex conjugation. 
Since = — 1 for spin 1/2 electrons, it follows from 
Kramers' theorem that every Bloch state at wavevector 
k is degenerate with a time reversed Bloch state. There- 
fore, the energy bands come in pairs, which are degener- 
ate at the two time reversal invariant points k* = and 
TT, as shown in Fig. 2. Note that in the presence of spin 
orbit interactions these bands can not be labeled with 
spin quantum numbers. 

In section IIIA we related the charge polarization as 
the sum of the Wannier centers of all of the bands. 
Kramer's theorem guarantees, however that the Wan- 
nier states come in Kramer's degenerate pairs, in which 
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each pair has the same center. The idea is therefore to 
keep track of the center of one of the degenerate Wannicr 
states per pair by defining a "partial polarization" . This 
will contain more information than Eq. H3.4|l . which is 
the sum over both states. 

For simplicity we assume that there are no degenera- 
cies other than those required by time reversal symme- 
try. Therefore, the 2N eigenstates may be divided into 
N pairs which satisfy 



,11 



-k,al 



(3.10) 



where a = l,....iV. The second equation follows from 
the first, along with the property 0^ = — 1. As shown 
in Fig. 2 these bands are defined continuously at the 
degeneracy points k* — 0, tt. This representation is not 
invariant under the general U{2N) transformation \6.'6\i . 
However, that invariance will be restored below. 

We define Wannier states associated with these two 
sets of bands along with the corresponding Wannier cen- 
ters. By analogy with (|3.4|) the partial polarization as- 
sociated with one of the categories s — I or II may then 
be written 



where 



A'{k) 



«a|VfeK„). 



(3.11) 



(3.12) 



The partial polarizations are clearly invariant (up to a 
lattice translation) under changes in the phases of luj. ^) 
and \ul^^). However, they appear to depend on the ar- 
bitrary choice of the labels / and // assigned to each 
band. We now show that the partial polarizations (|3.11(l 
can be written in a form that is invariant under a gen- 
eral U{2N) transformation of the form H3.3|l . To make 
this invariance explicit for Pj we treat the portions of the 
integral for positive and negative k separately. 



P^ = -^ / dk \A\k) + A\-k) 
27r Jo 



(3.13) 



For the second term we use the time reversal constraint 



(jXTni) along with the fact that (eu(^„|Vfe|e 



,11 



,11 



,11 



to write 



gauge invariance. The second term, can be rewritten in 
a suggestive manner by introducing the U{2N) matrix 
which relates the time reversed wavefunctions. 



(3.16) 



In the representation (|3.1UI) Wmn is a direct product of 
two by two matrices with e'-^'='° and —e'^^-*'-" on the off 
diagonal. At fc = and k — tt Wmn is antisymmetric. 
An antisymmetric matrix may be characterized by its 
Pfaffian, whose square is equal to the determinant. We 
then find that 



Pi[w{7T)] 

Pf[w(0)] 



exp[jy](X7r,a - XO,a)]- 



(3.17) 



Thus, the second term in H3.15|l can be expressed in terms 
of Pf[w]. This leads to 



P' = — 
2-K 



dkA{k) + Hog 



Pf[w(7r)] 
Pf[w(0)] 



(3.18) 



Using the identity Pf[XAX^] = Det[A:]Pf[A] it can 
be shown that under the U{2N) transformation l|3.3|l 
Pf[w] Pf[w]Det[C/]. Both terms in (|TTH|) are thus 
clearly SU{2N) invariant. Moreover, under a U{1) trans- 
formation the two terms compensate one another, so P^ 
is U{2N) invariant. Like the charge polarization (|3.4|l . 
P^ is only defined modulo a lattice vector. This is re- 
flected in the ambiguity of the imaginary part of the log 
in (|3.18l) as well as the dependence of gauge transfor- 
mations where the phase of \uk.n) advances by 27r for 

< /C < TT. 

A similar calculation can be performed for P^^ , and 
it is clear from time reversal symmetry that P^^ = P^ 
modulo an integer, reflecting the Kramer's pairing of the 
Wannier states. From H3.4|) and (|3.11|) the charge polar- 
ization is given by the sum of the two partial polariza- 
tions: 



Pp^P' 



P 



II 



(3.19) 



We now define the time reversal polarization as the dif- 
ference: 



Pe =P'-P 



II 



= 2P^ 



(3.20) 



A'i~k) = A"{k) - VkXk,a- (3.14) This then has the form 

It then follows that 



Pb 



P' 



1 

2^ 



/ dkA{k) ~T{x.., 
Jo 



Xa.a) 



1 

2^ 



dkA 



dkA + 2ilog ' ^ ^ " 



(3.15) 



(3.21) 

This may be written more compactly in terms of w„m as 



The first term is expressed in terms of the Berry's connec- 
tion A = A^ + AJ^ ■ However, since the path of integra- 
tion is not closed, the second term is necessary to preserve 



1 

2jxi 



[ dkTr[w^Wkw] - 21og 
Jo 



Pf[w(7r)]' 

pfRo)!, 



(3.22) 
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This can further be simphfied by noting that the first 
term gives the winding of the U{1) phase of between 
and TT. Thus, 



1 

2Td 



dfcVfelogDet[w(fc)] - 21og 



Pf[w(7r)] 



Pf[w(0)] 

(3.23) 

Since Det[i(j] = Pf[i(;]^ this quantity is an integer, and 
the due to the ambiguity of the log, this integer is only 
defined modulo 2. Even and odd integers are distinct, 
however, and determine whether Pf is on the same 

branch or opposite branch of •\/Det[w(fc)] at A: = and 
TT. An alternative way of writing it is thus. 



(-1) 



Pe _ 



^Det[w(0)] v/Det[w(7r)] 
Pf[w(0)] Pf[u>(7r)] ' 



(3.24) 



where the branches of ±-\/Det [w] are chosen such that 
the branch chosen at fc = evolves continuously along 
the path of integration in (|3.23|) into the branch chosen 
at fc = TT - eliminating the ambiguity of the square root. 

Eqs. H3.21I3.24|I are among the principle results of this 
paper, and can be regarded as a generalization account- 
ing for time reversal symmetry of the Berry's phase for- 
mulation of the charge polarization?-. The Z2 time re- 
versal polarization Pg defines two distinct polarization 
states. In the next section we will argue that the value of 
Pg is related to the presence or the absence of a Kramers 
degenerate state at the end of a finite system. As is the 
case for Pp, the value of Pg is not meaningful by itself, be- 
cause a gauge transformation e^'^lw^) changes its 
value. Equivalently, the presence or absence of a Kramers 
degeneracy at the end can not be determined from the 
state in the bulk, since it will depend on how the crys- 
tal is terminated. Nonetheless, the two values of Pg are 
topologically distinct in the sense that the value of Pg can 
not be altered by a continuous change in the Hamiltonian 
which preserves time reversal symmetry. However, in the 
next section we will argue that an adiabatic change in the 
Hamiltonian which preserves time reversal symmetry at 
the end points - but not in between - leads to a well 
defined change in Pg. This change defines a topological 
classification of distinct pumping procedures. 



C. Z2 Invariant 

In the previous subsection we focused on a time rever- 
sal invariant Hamiltonian, which occurs at t = and T/2 
in our pumping cycle. We now consider the continuous 
evolution of the Hamiltonian through the cycle and show 
that the change in the time reversal polarization which 
occurs in half the cycle defines a Z2 topological invariant 
which distinguishes a Z2 spin pump from a trivial cycle. 

This physical meaning of this invariant is easiest to see 
pictorially by considering the shift in the Wannier centers 
in the course of one cycle. Fig. 3(a) depicts the centers 
of the occupied Wannier orbitals as a function of t. At 




FIG. 3; (a) Schematic diagram showing the evolution of the 
centers of the time reversed pairs of Wannier states as a func- 
tion of t. Between t — and t = T/2 the Wannier states 
"switch partners", resulting in the appearance of unpaired 
Wannier states at the end. (b) The torus defined by k and t, 
with the four time reversal invariant points Fi connected by 
paths C12 and C34. 



t = 0,T/2 and T, time reversal symmetry requires that 
the Wannier states come in time reversed pairs. How- 
ever, in going from t = 0tot = T/2 the Wannier states 
"switch partners" . In this process, the time reversal po- 
larization, which tracks the difference between the posi- 
tions of the time reversed Wannier states changes by one. 
In addition, this switching results in the appearance of 
an unpaired occupied Wannier state at each end. Since 
the Wannier states come in pairs, there must be twofold 
Kramers degeneracy associated with each end - resulting 
in a total degeneracy of four. 

When the system evolves from t = T/2 to t = T, 
there is another switch, and the time reversal polar- 
ization returns to its original value. However, since 
H[t] — QH[T — t]0~^, the system with open ends does 
not return to its original state at t = but its ends are in 
an excited state because of the level crossing at t = T/2. 

We now relate the occurrence of this non trivial pump- 
ing cycle to a topological property of the bulk ground 
state as a function of t. We thus consider the change in 
the time reversal polarization, Pg(t) between t = and 
t = T/2. Note that though Pg is not gauge invariant, the 
change in Pg is gauge invariant. This difference 



A = Pg{T/2) ^ Pg{0) mod 2 



(3.25) 



defines a Z2 topological invariant which characterizes the 
mapping from the torus defined by fc and t to the wave- 
functions: \uk,nit))- From H3.24|l we may write this in- 
variant as 



v/Dct[^;(r,)] 

pf[t.(r,)] ■ 



(3.26) 



Here F^ are the four "time reversal invariant points" on 
the torus shown in Fig. 3(b). The branches of the square 
root are chosen as in (|3.24l) by continuously evolving 
■\/Det[w(fc, t)] along the paths C12 and C34. In order to 
apply this formula, it is crucial for the wavefunctions to 



8 



be defined continuously on the torus. It is always possible 
to find such smoothly defined wavefunctions via a trans- 
formation of the form because the Chern number - 
which is the obstruction to doing so - is equal to zero. 

In the Appendix we will relate different mathematical 
formulations of this invariant. We will first show that 
it can be interpreted as an obstruction to defining con- 
tinuous wavefunctions provided an additional constraint 
relating the wavefunctions at time reversed points is en- 
forced. This leads to a different formula for the invariant, 
which can be expressed in terms of the Berry's curvature 
^ and the Berry's connection A. We will then prove 
that H;-{.26(l is equivalent to the formula for the invariant 
introduced in Ref. 



IV. ELECTRON INTERACTIONS AND 
BOSONIZATION 



The preceding discussion has focused on non interact- 
ing electrons. An important question is therefore whether 
or not these ideas apply to interacting systems. The pres- 
ence or absence of a groundstate Kramers degeneracy as- 
sociated with the ends of a finite interacting time reversal 
invariant system is clearly a well posed yes or no ques- 
tion. This suggests that the time reversal polarization is 
a well defined quantity - at least for non fractionalized 
phases for which the groundstate with periodic bound- 
ary conditions is non degenerate. Therefore, we believe 
the topological distinction of the Z2 pump is still present 
with interactions. 

Calculating the time reversal polarization for interact- 
ing electrons is more subtle than for non interacting elec- 
trons. One possible approach would involve character- 
izing the entanglement entropy, as in Ref. ISlL which is 
sensitive to the presence or absence of end states. In this 
section we adopt a simpler approach by studying an inter- 
acting version of the model introduced in section II using 
abelian bosonization. We find that bosonization provides 
a natural description of the time reversal polarization. 

We begin with a continuum version of l|2.1|) described 
by the Hamiltonian density, 

7? = (ivpT'd^ + h^tT^'a' + AtstT^ + icso ■ ar') ip. 

(4.1) 

Here if^aa is a four component field, where the left and 
right moving fields a = L, R are specified by the eigenval- 
ues of r^j and the spin a —H by a^^. We now bosonize 
according to 



with 



1 



\J2-KXc 



(4.2) 



where short distance cutoff. Define charge/spin 

variables so that 01 /^ct = 0pa i ^aa-, and charge/current 
variables (with ^ — p,a) as 4>^l/r — ^fj. ^ S^. These 
obey [dx6^,{x), ipf,' {x')] = i{n/2)6^^'6{x - x'). 

The bosonized Hamiltonian then has the form (|2.1|) 



Vh = 7/^sin26'pSin26icr, 

ilTXr 



Vt 



Atst 

2-KX, 



■ sin26'p cos 26*0-, 



(4.3) 
(4.4) 

(4.5) 



and 



Vsc 



TT 2,t:Xc 



sin 29 cr COS2tfa~ 



2'KX, 



■ sin 29 „ sin 2Lp„. 
(4.6) 



In the absence of the spin orbit term, the spin sector of 
this Hamiltonian (when Op is pinned at 7r/4) is equiva- 
lent Shindou's model^. This Hamiltonian describes an 
insulating phase in which both Op and 0^ are pinned. 

First focus on the case hst = where the Hamiltonian 
is time reversal invariant. If we choose a gauge such that 
Q'ij)<d~^ = T^a-y^'*, the behavior of these operators under 
time reversal can be deduced: 



eOpS'^ = 



'pi 



OippQ ^ = -ipp, 

e^^e-^ = if„ + 7T/2. (4.7) 



The time reversal invariance of the Hamiltonian when 
hst = can easily be verified. It is now straightforward to 
consider time reversal invariant interaction terms, such as 
(dxOp)'^, (dxOa)'^, cos40p, cos46'ct, cosiipa, etc. Provided 
these interaction terms (as well Vso defined above) are 
not too large, the system will retain its bulk gap and be 
in a phase in which Op and 6a- are pinned. 

We now identify the time reversal polarization with 



Pg = 2da/TT mod 2. 



(4.8) 



The apparent dependence of Pg on the spin quantization 
axis is an artifact of abelian bosonization. In fact, Pg is 
SU{2) invariant. This can be seen by noting that global 
spin rotations are generated by 5^ ~ / dxdxOa and ~ 
/ da;exp±2i0CT- The latter obeys [6*^,5'^] = ztTrS*^, so 
that [Pe, S^] = 0. It can further be seen that even in the 
presence of spin nonconserving terms in Vso as well as 
the interaction terms discussed above, [Pg, H] — 0. Since 
QPgQ~^ = —Pe mod 2 there are two distinct possible 
values for the time reversal polarization: (Pg) = or 1. 
Thus Pg can be used to classify time reversal invariant 
insulating states. 

Consider a finite system with ends. We now argue 
that the value of P^ determines the presence or absence 
of Kramers degenerate states at the ends. The end of a 
one dimensional system at x = must be characterized 
by a boundary condition for 6cr{x ~ 0). Time reversal 
symmetry limits the possible values to Oa{x = 0) = mr. 
The value of n, however, depends on how the lattice is 
terminated. First suppose that n = 0. Then, when P^ = 
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1, the pinning of 9a in the bulk is not consistent with 
the boundary condition. The closest it can be is {6a) — 
±7r. Thus, near the end there must be a kink of ±7r in 
9a at the end. Time reversal symmetry requires these 
two possibilities to be degenerate, so there is a Kramers 
degeneracy of two at the end. On the other hand, when 
Pa = 0, the bulk energy gap is "consistent" with the 
boundary condition, allowing for 9a{x) = everywhere. 
The groundstate in this case is unique. 

We thus conclude that bosonization provides an alter- 
native approach for formulating the time reversal polar- 
ization in terms of 9a , just as it allows for a formulation 
of the charge polarization Pp = 9p/'K. This suggests that 
the topological distinction of the Z2 spin pump remains 
in the presence of electron interactions. 



V. DISCUSSION 



A. Relation to Quantum Spin Hall Effect 




Gap 



0(1/R)j_ 

T 







1/2 




FIG. 4: Schematic plots of the many particle eigenstates of a 
cylindrical quantum spin Hall system of radius i? as a func- 
tion of the magnetic flux threading the cylinder. In each case 
there is a Kramers degeneracy when the flux $ is equal to 
h/2e. (a) There exist gapless edge excitations, whose energy 
level spacing goes to zero for R — > 00. (b) There is an edge ex- 
citation gap which remains finite for R 00, but the ground 
state is doubly degenerate. 



The quantum spin Hall phase introduced in Ref. ^3 is 
a phase of a two dimensional electron system. In a man- 
ner analogous to Laughlin's construction for the quantum 
Hall effect'^?, this phase, when compactificd onto a cylin- 
der, defines a Z2 pump of the sort studied in this paper. 
In this section we outline the implications of the present 
work for the quantum spin Hall effect. We begin by re- 
lating the Z2 index introduced in section III to the index 
which distinguishes the quantum spin Hall phase from a 
band insulator. We then discuss the presence or absence 
of gapless edge states in the quantum spin Hall effect. 
Finally, we comment on an alternative topological char- 
acterization of the quantum spin Hall effect in terms of a 
"Chern number matrix" that has recently been proposed 
by Sheng et ali^. 

1. Z2 classification of quantum spin Hall phase 

For non interacting electrons, the electronic phase of 
a two dimensional system with a bulk gap is character- 
ized by the wavefunctions defined on the Brillouin zone 
torus, \un{kx, ky)). The relationship between the one di- 
mensional Z2 pump and the two dimensional quantum 
spin Hall effect can be established by the identification 
of (fc,t) with {kx,ky). Eq. H1.2|l then reflects the time 
reversal invariance of the two dimensional Hamiltonian. 
As we prove in the Appendix, the Z2 topological index 
introduced in Ref. 13 is equivalent to the Z2 index char- 
acterizing the pump. The considerations of this paper 
provide a natural physical interpretation of this index in 
terms of the change in the time reversal polarization in 
half of the cycle. 

In addition, our observation that the time reversal 
polarization is related to the presence or absence of a 
Kramers degeneracy at the end suggests that the Z2 
classification of time reversal invariant two dimensional 



ground states transcends the non interacting model for 
which it was derived. This means that the quantum spin 
Hall effect describes a phase that is distinct from a band 
insulator even in the presence of electron electron inter- 
actions. 



2. Edge states or not? 

In the regular quantum Hall effect, the topological 
structure of the bulk state guarantees the existence of 
gapless edge excitations. The non trivial Z2 invariant, 
however, does not guarantee gapless edge states. As 
shown by Wu et al.-^ and Xu et aliiS:, when the inter- 
actions at the edge are sufficiently strong the edge can 
undergo a transition which opens a gap. The consider- 
ations of this paper allow us to prove that either there 
are gapless edge excitations or there is a ground state 
degeneracy at the edge associated with the breaking of 
time reversal symmetry. 

To establish this proposition, consider the many par- 
ticle eigenstates of the quantum spin Hall phase on a 
cylinder as a function of the magnetic flux through the 
cylinder. When the radius R of the cylinder large, then 
the 0{1/R) change in the energy of the many particle 
eigenstates when the one half a flux is inserted will be 
much less than any energy gap. But in the quantum spin 
Hall state, the non trivial Z2 index requires the ground 
state to have a Kramers degeneracy at either $ = or 
<& — 4>q/2, but not both. Thus there are two possibilities 
as schematically illustrated in Fig. 4. Either there are 
edge states with energy 0{1/R) which become gapless 
for i? — > 00 or the ground state is degenerate for R 00 
and split by at most 0{1/R) by the magnetic flux. 

This required ground state degeneracy distinguishes 
the quantum spin Hall phase from that of a band in- 
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sulator. Unlike a band insulator, the quantum spin Hall 
state in a system with edges can not have a unique ground 
state with a gap for all excitations. 

3. Other proposed classifications of the quantum spin Hall 

effect 

We now comment on a different topological classifica- 
tion of the quantum spin Hall effect proposed by Sheng 
et al^ These authors defined a matrix of Chern num- 
bers by considering a system with a generalized class of 
periodic boundary conditions. Specifically, they consid- 
ered boundary conditions of the form $(..., -1- Lj , ... ) = 
exp(i6'")<i>(..., r^Q, ...), where Lj=a:,y define the periodic- 
ity and are independent phase twists for the up 
and down spins. They then characterized the topological 
classes of the groundstate wavefunction as a function of 
these phase twists, and defined a matrix of Chern num- 
bers, C"''^. 

Sheng et al. argued that this classification contains 
more information than the Z2 classification because it 
distinguishes quantum spin Hall states which belong to 
the same Z2 class. This can be illustrated by looking at 
the continuum version of the graphene model introduced 
in Ref. 0, described by the Hamiltonian 

H = [-iVF{(TxTzdx + O-ydy) + AsoCTzT^Sz] ip. (5.1) 

Here, in the notation of Ref. ^j, az describes the sub- 
lattice of the honeycomb lattice, Tz describes the two 
inequivalent valleys at the corners of the Brillouin zone 
and Sz describes the spin. When A^o is nonzero the sys- 
tem is in a quantum spin Hall phase and belongs to the 
nontrivial Z2 class. Sheng et ali^ argued that the sign of 
Aso defines two distinct phases which are distinguished 
by the matrix of Chern numbers. 

When Sz is conserved this is certainly correct, and 
the Chern number matrix can be viewed as independent 
Chern numbers for the up and down spins. However, 
when spatial symmetries are relaxed and spin is not con- 
served this distinction is no longer meaningful. The two 
phases discussed above are in fact the same phase because 
they can be continuously transformed into one another 
without closing the gap. Specifically, consider the more 
general spin orbit interaction which preserves the energy 
gap: 

azTzSz ^ (TzTz{s- n). (5.2) 

When the unit vector h is continuously varied from +2 
to —z the two "phases" are connected. Of course, the 
process of connecting these phases requires the breaking 
of the C3 lattice symmetry of graphene. But in general, 
disorder will break all spatial symmetries, so one can not 
rely on a spatial symmetry to protect a topological prop- 
erty. 

This presents a conundrum because the Chern ma- 
trix formulation distinguishes the two states with dis- 
tinct topological integers - even when the C3 symmetry 



is explicitly violated. What happens to these integers 
when the continuous path in Eq. I|5.2|l is adiabatically 
followed? The answer is that somewhere along the path 
the energy gap must vanish at the edge where the twisted 
spin boundary condition is imposed^'^. 

The spin phase twist imposed by Sheng et al. can be 
decomposed into a C/(l) part 6p = 9^ + 9^ and a "spin" 
part ^o- = 9^ — Oi- The spin phase twist 9^ is fundamen- 
tally different from 9p when the bulk Hamiltonian does 
not commute with the Sz ■ The boundary where the spin 
phase twist is imposed is physically different from the rest 
of the system, and the spectrum of the Hamiltonian will 
in general be different for different values of 9a- ■ In con- 
trast, the location of the charge phase twist 9p introduced 
by Niu and Thouless^, can be moved around by per- 
forming a local gauge transformation without changing 
the spectrum. Since the vector potential can be spread 
out over the circumference 2ttR of the torus the change 
in the spectrum due changing 9p can be at most 0{1/R). 
In contrast, if the spin phase twist is spread out over the 
circumference, changing 9a- changes the Hamiltonian by 
an amount of order 1. The spectrum need not be close 
to that of the physical Hamiltonian. 

We conclude that the additional topological structure 
implied by the Chern number matrix is a property of the 
boundary where the twisted phase condition is imposed 
rather than a property of the bulk two dimensional phase. 
The bulk quantum spin Hall effect is classified by the Z2 
invariant alone. 



B. Can the Z2 spin pump pump spin? 

Is the Z2 pump we introduced a spin pump? Since an 
isolated Z2 pump returns to its original state after two 
cycles, the simple answer to this question is no. How- 
ever, any functioning pump must be connected to reser- 
voirs into which the pump can pump. In this section we 
briefly consider the effect of connecting the Z2 spin pump 
to reservoirs. We conclude that the Z2 pump does pump 
spin, though the spin pumped per cycle is not quantized. 
Moreover, we argue that when the coupling to the reser- 
voirs is weak the Z2 topological structure of the pump is 
essential for a nonzero spin to be pumped. For stronger 
coupling, however, the Z2 structure is not essential. 

We consider a simple case where the reservoirs can 
be described by non interacting electrons with vanishing 
spin orbit interaction. We first suppose the coupling to 
the reservoir is weak, so that the level width T induced 
in the pump due the coupling is small compared to the 
energy gap A. However, we require the coupling T to 
be large compared to the pumping rate, as well as any 
inelastic scattering rate for the end states. In the limit 

VT, n/r^ < r < A, (5.3) 

the eigenstates of the pump maintain their integrity, 
though coupling to the reservoirs allows transitions be- 
tween different states. 
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FIG. 5: Evolution of one end of a pump that is weakly 
connected to a lead as shown in inset. For the excited 

state of the pump can relax by creating a odd under time 
reversal excitation in the lead, which generically will change 
the spin of the reservoir. 



for it could be a Kondo resonance where the Kramers de- 
generate end state becomes entangled with the reservoir 
electrons. 

It should be emphasized that the spin added to the 
reservoir is not a property of the bulk Hamiltonian of 
the pump, but rather it depends on how the pump is 
connected to the reservoir. The spin transferred to the 
reservoirs at the two ends of the pump need not be re- 
lated. Thus, one can not view the spin as being pumped 
along the length of the pump. However, the presence of 
the end state resonance, which follows from the change in 
the time reversal polarization, is a property of the bulk 
insulating state. In this sense, the Z2 pump is a pump 
for spin. 



As illustrated in Fig. 5, there is a point in every cycle 
t = (n + l/2)r, where the groundstatc of the pump be- 
comes degenerate. This degeneracy is due to the end 
states, which are in proximity to the reservoir. For 
t {n -\r 1/2)T the pump is in an excited state. Cou- 
pling to the reservoir, however, allows the pump to relax 
back to its groundstate. This relaxation, however, must 
involve a process in the reservoir which is odd under time 
reversal. Generically this will involve changing the spin 
of the reservoir. The spin added to the reservoir need 
not be quantized. It could even be equal to zero, but 
generically, it will be of order h. 

In this weak coupling limit it is clear that the Z2 struc- 
ture of the pumping cycle is essential because it guaran- 
tees the level crossing in the end states. If the end states 
did not cross, then there would be no transitions, and the 
spin in the reservoirs would be unchanged after a com- 
plete cycle. However, finite coupling between to leads 
relaxes this requirement. Suppose that the time reversal 
symmetry is weakly broken at t = T/2, so that there is 
a small anticrossing of magnitude 5. In this case the Z2 
character of the cycle is lost. But if i5 <C F, then the 
states have no way of "knowing" about the anticrossing, 
and the pump proceeds as if (5 = 0. 

This reflects the fact that spin can be introduced into 
a reservoir which is connected to an insulating material 
when the insulator is deformed through a periodic cy- 
cle. The spin injected can be expressed in terms of the 
unitary reflection matrix f(t) for electrons at the Fermi 
energy in the reservoir—, which in general depends on 
the Hamiltonian H{t) of the insulator, 

AS= — ((dt Trf^ft^l. (5.4) 
2ni J ^ dV ^ ' 

In general, this quantity is non zero. The difficulty is 
coming up with a cycle in an insulating material for which 
is not very small. The Z2 pump accomplishes this by 
guaranteeing that there is a resonance in the reflection 
matrix, which occurs when the Kramers degenerate end 
state appears. Note that this resonance need not involve 
charge transfer between the reservoir and the insulator. 
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APPENDIX A: EQUIVALENT FORMULATIONS 
OF THE Z2 INVARIANT 

In this appendix we relate different mathematical for- 
mulations of the Z2 invariant A. Our starting point is 
Eq. along with Eqs. p. 2 113. 24(1 which express the 

invariant in terms of the change in the time reversal po- 
larization between t = Q and t = T/2. We will first show 
that A can be interpreted as an obstruction to defin- 
ing wavefunctions continuously, provided a time reversal 
constraint is enforced. This will lead to a formula for A 
in terms of the Berry curvature and Berry connection, 
which will be shown to be equivalent to (|3.25|) . We will 
then show that (|3.25|l is equivalent to the Z2 invariant 
proposed for the quantum spin Hall effect in Ref. IT^ 

We will use a notation appropriate for the Z2 pump- 
ing problem and consider Bloch wavefunctions defined 
continuously on the torus defined by — tt < k < tt and 
< t < T. For the two dimensional quantum spin Hall 
effect wc should identify with k and ky with 2TTt/T. 



1. Z2 invariant as an obstruction 

It is well known that a nonzero value of the Chern 
invariant is an obstruction to smoothly defining the 
wavefunction throughout the entire torua^^*^. Instead, 
wavefunctions must be defined on overlapping "patches" 
which are related to each other by a gauge transforma- 
tion called a "transition function" . The Chern number is 
then related to the winding number of the phase of the 
transition function around a non contractable path. For 
the problem studied in this paper, the Chern number is 
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zero, so there is no obstruction to finding a transforma- 
tion of the form which makes the wave functions 
smoothly defined on a single patch. However, we will 
show in this section that if we enforce the time reversal 
constraint 



\uii-k,-t)) = e\ui\k,t)) 

\ui'i-k,-t))=-Q\uiik,t)), 



(Al) 



then a nonzero value of the Z2 invariant is an obstruction 
in a manner precisely analogous to the Chern number. 
This constraint means that the gauges for the wavefunc- 
tions at ±(k,t) are not independent. At the four time 
reversal invariant points (fc,t) = F^, the allowed trans- 
formations of the form are restricted to be symplec- 
tic, Umni^i) G Sp{N). That a nonzero value of the Z2 
invariant A is inconsistent with this constraint is easy to 
see because it implies that Det[w{k, t)] = 1 for all k and 
t and Pf[w(ri)] = 1, so Eq. trivially gives A = 0. 

We will now relate the Z2 invariant to the winding of 
the phase of transition functions relating the wavefunc- 
tions on different patches. In addition to establishing the 
connection between the Z2 invariant and the Chern in- 
variant, this approach will derive a formula for the Z2 
invariant which expresses it in terms of the Berry's con- 
nection and curvature. The similarity between the Z2 
invariant and the Chern invariant has been emphasized 
by Haldane^. The formulation of the Z2 invariant as 
an obstruction has also been discussed by Roy'^^, though 
that work did not establish a formula for the invariant. 

Suppose that we have wavefunctions obeying I)A1|I de- 
fined smoothly on two patches in the torus labeled A and 
B in Fig. 6. In patch A the wavefunctions \u^{k, t))A are 
smoothly defined everywhere in the upper left and lower 
right quadrants of Fig. 6, while for patch B \u%{k,t))B 
are defined in the upper right and lower left quadrants. 
In the overlapping regions these different wavefunctions 
are related by a U{2N) transition matrix 



.{k,t)). 



^rnn I "i: 



(fc,t))j 



(A2) 



where m and n are shorthand of s and a. Consider the 
change in the C/(l) phase of t^^ around the closed loop 
9ti in Fig. 6. 



27ri Jar, 



(A3) 



This will clearly be an integer because it is equal to the 
winding number of the phase of Det[t'^^] around the loop 
dri. If D is non zero and can not be eliminated by a 
gauge transformation, then there is an obstruction to 
smoothly defining the wavefunctions on a single patch. 
In what follows, we show that D mod 2 is precisely equal 
to the Z2 invariant defined in this paper. 
From ljA3|l . we may write 



D 



1 

2^ 



(A4) 



T/2 



A 



'1/2. 



B 



B 



A 



-71 







71 



FIG. 6: The torus defined by k and t divided into two patches 
A and B. The boundaries of the regions ri and T2 are shown 
as well as the boundary of the shaded region, ri/2. 



where A'^ = X^nA^^"!^!"")^ ^^"^ likewise for A^ . Since 
\un)A is smoothly defined in the interior of ri we may 
write it in terms of the Berry's flux, 



(A5) 



This can not be done for |u„) b which is not necessarily 
defined continuously inside t\. However, it can be related 
to the Berry's fiux through T2, 

/ d£-A^^~(f dl-A^+(f di-A^ (A6) 

J dri <J dr2 '^Q'^'i/i 

= ~ j drT^ + (f dl-A^. (A7) 

Combining these we thus find the winding number for 
the transition function can be expressed as an integral 
involving the Berry connection and the Berry curvature. 



D 



1 

2^ 



dt-A- / drT 



mod 2. 



(A8) 



9ri 



The patch labels can be safely removed because is 
gauge invariant, and the line integral is gauge invariant 
modulo 2. It is essential that the time reversal constraint 
(|A1(I be enforced for this equation to have meaning. If 
not, then a gauge transformation on patch B can change 
the line integral by 1, making the formula vacuous. When 
(|A1I) is enforced, an odd value of D can not be gauged 
away because the phases of \u^) and can not be in- 
dependently changed. Thus, D — 1 mod 2 presents an 
obstruction to defining wavefunctions on a single patch. 

We now show that this winding number is precisely 
the same as the invariant A. To this end, we rewrite A 
in terms of the partial polarization defined in 13.111 



13 



Using Pg = 2P^ — Pp we have 

A = 2 (P^(T/2) - P'{0)) - (Pp(r/2) - PpiO)) mod 2. 

(A9) 

When (|Aip is enforced Eq. (|3.18|) for the partial polar- 
ization implies that 



2(P^(r/2) -P^(0)) = 

Eq. H3.4|l shows that 

Pp{T/2) - Pp{0) 



27r 



di ■ A. (AlO) 



9ti, 



27r 



drT. 



(All) 



Combining the two terms thus establishes that A — D. 
The two terms in (|A8(I thus acquire physical meaning: 
The line integral gives twice the change in the partial 
polarization between t = and t = r/2, while the surface 
integral gives the change in the total polarization. 



2. Zeros of the PfafHan 

In Ref. 13, the Z2 invariant was introduced by consid- 
ering the matrix elements of the time reversal operator, 



mij{k,t) = {ui{k,t)\Q\uj{k,t)). 



(A12) 



This should be contrasted with the matrix Wij{k,t) in- 
troduced in section III, which can be generalized as a 
function of t to be 



Wij{k,t) = {ui{~k,~t)\Q\uj{k,t)). 



(A13) 



At the four time reversal invariant points (fc, t) = 
ri.2,3,4 = (0, 0), {tt, 0), (0, r/2), {tt, T/2), w,j and m,, co- 
incide, but in general they are different. Wij is uni- 
tary with |Det[w]| = 1, while rriij is not unitary. Since 
= —1, rriij is antisymmetric. The Pfafhan of m is 
therefore defined for all k and t. In Ref. ^3 we argued 
that the Z2 invariant could be determined by counting 
the number of zeros of the Pfafhan in half the torus, mod- 
ulo 2. 

To establish the equivalence of this with (|3.25(l . we 
begin by rewriting the time reversal polarization Pg in 



terms of Pf[TO(fc,t)]. The key observation to be made is 
that 



Det[w(fc,t)] 



Pf[TO(fc,t)] 

Vi[m{-k,-t)\* 



(A14) 



which can be proved by noting that m{~k^~'t) — 
w{k,t)m{k,t)*w{k,t)'^ and using the iden- 
tity Pf[XAX^] = Det[A:]Pf[A]. Introducing 
p{k) = Pf[m{k,t*)] for t* = 0,T/2 it follows that 
logDet[w {k,t*)] = log p{k) + log p{—k). Thus we may 
rewrite (j3.23|l as 



Pa = 



1 

27ri 



dfcVfc[log p{k) -f log p{-k)] - 21og 



p{tt) 



MO) 
(A15) 

where we have used the coincidence of w and m at k = 
and TT. This may be simplified further by changing 
variables A; — > — fc in the middle term and writing the 
last term as an integral from to tt. This gives 

1 r 

Pe^T^ I dfcVfelogPf[m(fc,0)] mod 2, (AI6) 



where the integral is now over the closed loop t = Q, 
— TT < A; < TT. This expression is only defined modulo 2 
because of the ambiguity of the imaginary part of the log 
in jSUl). 

Thus, we have established that Pg is given by the phase 
winding of the Pfaffian, p(k) around the ID Brillouin 
zone modulo 2. While this quantity is not gauge invari- 
ant, the change in it due to continuous evolution between 
t* = and t* = r/2 is gauge invariant. This defines the 
Z2 topological invariant, which, as in Ref. Il3l may be 
written 

A = — / dt- VlogPf [m(/c, t)] mod 2, (A17) 
27^* /an/, 

where dri is the boundary of half the torus defined by 
-TT < fc < TT and < t < T/2 (see Fig. 6). If Vi\m(k,t)\ 
has point zeros, then this quantity counts the number of 
zeros in rx/2 modulo 2. 
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